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Abstract: We calculate the amplitude of gg — * bb production for the 
colour singlet, J z = di-gluon state at 0(a|) order. We consider the cancel- 
lation, and a realistic cutoff, of the infrared divergent terms. We show that 
the one-loop radiative QCD contributions effectively reduce the Born level 
result for the central exclusive bb cross section at the LHC. This process is 
essentially the only irreducible QCD background to the exclusive H — > bb 
signal. 
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1 Introduction 



The search for, and identification of, the Higgs boson(s) is one of the main 
goals of the LHC. Once the Higgs boson is discovered, it will be of primary 
interest to determine its spin and parity, and to measure precisely its mass 
and couplings, in particular the Hbb Yukawa coupling. An important con- 
tribution to the comprehensive study of the Higgs sector can be provided by 
central exclusive diffractive (CED) Higgs boson production, 

pp ^ p® H ®p, (1) 

where the © signs denote the presence of the large rapidity gaps (LRG); see, 
for example [!]-[!]■ The process is sketched in Fig. [TJ In such an exclusive 
process there is no hadronic activity between the outgoing protons and the 
decay products of the central (Higgs) system. The bb decay mode of the 
Higgs is especially attractive, since its observation would allow a detailed 
study of the MSSM Higgs sector [5]- [TJ. Indeed, for certain BSM scenarios, 
it may become the Higgs discovery channel [5j [8] . The experimental study of 
central exclusive Higgs production is one of the key theoretical motivations 
behind the FP420 project [HI EI] which proposes to complement the CMS 
and ATLAS experiments at the LHC by installing additional forward proton 
detectors 420 m away from the interaction region. 

To avoid the production of new secondaries across the LRG, the colour 
flow caused by the active gluons (which participate in the gg —>■ H fusion) 
must be screened by another, more soft, gluon, see Fig. [TJ Thus, in the 
exclusive process the central system (the Higgs boson) is actually generated 
by a specific gg pp gluon state, where the PP superscript is to indicate that 
that each hard (active) gluon comes from colour-singlet t-channel (Pomeron) 
exchange. Moreover, the presence of the second t-channel screening gluon 
leads to an additional loop integration. The integration over the transverse 
momentum in this loop results in a J z = 0, CP-even selection rule [TT[ [T2] . 
Here J z is the projection of the total angular momentum along the proton 
beam axis. 

The J z = selection rule is one of the major reasons why central exclu- 
sive production is so attractive for Higgs studies. First, it readily permits a 
clean determination of the quantum numbers of the observed new resonance, 
which will be dominantly produced in a scalar state. A second direct conse- 
quence of this rule is the unique possibility to study directly the Hbb Yukawa 
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Figure 1: The mechanism for central exclusive Higgs production, pp —> p © 
H (Bp. The Higgs is produced by the fusion of two so-called active gluons. The 
screening gluon on the left is required to ensure a colour neutral system across 
the rapidity gaps. We do not show the additional screening corrections which 
must be included to ensure the survival of the rapidity gaps from population from 
secondaries resulting from soft rescattering. 



coupling of the Higgs-like bosons in central exclusive processes. The poten- 
tially copious 5-jet (QCD) background is controlled by a combination of the 
J z = selection rule [TTJ [TJ] (which strongly suppresses the leading-order 
QCD bb production), the colour and spin factors and, finally, the excellent 
resolution of the missing mass to the measured outgoing protons^. It is the 
possibility to observe directly the dominant bb decay mode of the SM Higgs 
with Mh ^5 140 GeV that first attracted attention to exclusive production 
at the LHC As is well known, the direct determination of the Hbb coupling 
appears to be very difficult for other search channels at the LHC. 

We emphasize, that for forward going protons at the LHC, the Higgs 
signal is produced by gluons in a J z = state whereas the LO QCD back- 
grounds are primarily initiated by initial states with \J Z \ = 2. The J z = 
background contribution is suppressed for large angles by a factor ml/E^, 
where E? is the transverse energy of the b and b jets, see for example [T31 [T4] . 
As discussed in [15j H3] , the physical origin of this suppression is related to 
the symmetry properties of the Born helicity amplitudes M x q \ q describing 

1 Current studies [9] show that the missing-mass resolution a will be of the order 1% for 
a 140 GeV Higgs assuming that both protons are detected at 420 m from the interaction 
point. The equality of the accurate missing-mass measurement of Mh with its mass 
determined from its decay products allows the background to be considerably suppressed. 
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the binary background process 

9(^1, Pa) + 9(^2,Pb) -> q(\,Pi) + q(Xg,p 2 ). (2) 

Here, the Aj label the helicities of the incoming gluons, and X q and A^ are the 
(doubled) helicities of the produced quark and antiquark. The p's denote the 
particle four-momenta (p A = p 2 B = 0, p\ 2 = m 2 ), with Pa+Pb = Pi +P2 and 
s — (pa +Pb) 2 - It was shown in [13] that for a colour-singlet, J z = 0, initial 
state, (Ai = A2 = A) the Born quark-helicity-conserving (QHC) amplitude 
with Xq = —X q vanishes^ 

<r Aa = o- (3) 

For the quark-helicity-non-conserving (QHNC) amplitude for large angle pro- 
duction we have 

M A A ,f ' ~ O {^\ Mtt\ (4) 

where the amplitude on the right-hand-side displays the dominant helicity 
configuration of the LO background process. 

The main sources of background to the exclusive H —> bb production 
at the LHC were discussed in detail in (31 [HI El [18]. It was shown that 
all backgrounds are strongly suppressed and controllable and, in principle, 
can be further reduced by the appropriate optimized cuts on the final state 
particle configurations in such a way that the signal-to-background ratio S/ B 
is of order 1 (or may be even better for MSSM Higgs [3 El [7]). 

Among all the QCD backgrounds, the m 2 / .^-suppressed di-jet bb pro- 
duction is especially critical, since it is practically the only one irreducible 
background source which cannot be decreased, either by improving the hard- 
ware (as in the case, when the prolific gg pp — > gg subprocess mimics bb 
production, and the outgoing gluons are misidentified as b and b jets [3]) or, 
for example, by cuts on the three-jet event topology (as in the case of large- 
angle gluon radiation in the process gg(J z = 0) — > qqg, discussed in [T4]). 
Therefore, the precise calculation of the QHNC background contribution be- 
comes of primary importance for quantifying the signal-to-background ratio, 

2 It is worth noting that in the massless limit Eq. ([3]) holds for any colour state of 
initial gluons. This is a consequence of the general property, that the non-zero massless 
tree-level amplitudes should contain at least two positive or two negative helicity states, 
see for example [16]. It is an example of the more general Maximally- Helicity- Violating 
amplitude (MHV) rule, reviewed for example in [17] , 
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and the accurate evaluation of the statistical significance of the H — > bb sig- 
nal. In [3J E] it was estimated that the contribution from this source was 
B/S ~ 0.4, when the Born formula for the binary cross section was used. 
The same LO result for the background subprocess was used for the evalu- 
ation of a statistical significance of the MSSM Higgs boson signal in [6j [7J. 
However, it is known ([HUE]) that higher-order QCD effects may strongly 
affect the LO result. 

First, there is a reduction coming from the self-energy insertions in the 
6-quark propagator, that is from running of 6-quark mass from m^mfc) to 
its value, mb(M#) < m^mb), at the Higgs scale. Here mb(/x) is the running 
6-quark mass in the MS scheme [19]. Secondly, in our case, where m& <C Ms, 
double-logarithmic corrections of the form (pis/it) In 2 (Mh/thi,) are poten- 
tially very important. These are related to the so-called non-Sudakov form 
factor, F q , in the cross section which arises from the virtual diagrams with 
gluon exchange, see [13] GO] 121] B2j . 

For the case of the ■y~f(J z = 0) — > bb process the complete one-loop result 
was first calculated in [20]. For the photon-photon reaction, the double- 
logarithmic (DL) asymptotics for F q has the form 



F q {L m ) = EM^J (5) 



71 

with L m = ln(Mf//ra&), c = 1 and c\ = —8 [T3], [20] . The second (negative) 
term in (jSJ) is anomalously large, and dominates the Born term for M# > 100 
GeV. Clearly, this dominance undermines the results of any analysis based 
on the Born approximation^. The physical origin of this non-Sudakov form 
factor was elucidated in [13], where the explicit DL calculation at the two- 
loop level was performed. It was also shown that the two-loop calculation 
should be sufficient for a reliable evaluation of the DL effects. This was 
confirmed by a more comprehensive all-orders study [2~Tl [22] . 

Recall that in the photon-photon case, the two-loop expression for the F q 
is p3] 

F g (L m ) = (l-3^) 2 + ^fl + ^) , (6) 



C 



F 



with 



as. 



T — — C F L m , (7) 

7T 



3 For gg pp fusion, rather than 77 fusion, we might, a priori, anticipate that colour 
factors could make the coefficient c\ even larger. 
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where Cp = (N% — 1)/2N C , Ca = N c and N c = 3 is the number of colours. 
The non-Sudakov logarithms come from kinematical regions where one of 
the quark propagators is soft. As well known, there are other DL effects 
(the well known Sudakov logarithms [23]) which arise from virtual soft gluon 
exchange. As discussed in [13], in the case of quasi- two-jet configurations, the 
Sudakov and non-Sudakov effects can be factorized to good accuracy, since 
they correspond to very different virtualities of the internal quark and gluon 
lines. For final state radiation, the Sudakov effects can be implemented in 
parton shower Monte Carlo models in the standard way. For the gg pp initial 
state, the Sudakov factors are explicitly incorporated in the unintegrated 
gluon densities, see [U [2] • 

Unfortunately, from a phenomenological perspective, it is dangerous to 
rely on the DL results, since experience shows that formally subleading (SL) 
corrections may be numerically important. That is why an accurate evalua- 
tion of the QCD bb background to the central exclusive process 

pp->p®(H->bb)®p (8) 

requires, first of all, the calculation of the exact one-loop correction to the 
Born gg pp — > bb amplitude, that is to the hard subprocess gg — > bb in a 
colour-singlet J z = initial state. 

The calculation is described in Section 2, where the result is presented 
as the sum of standard integrals corresponding to the box, triangle and self- 
energy Feynman diagrams with the scalar propagators. In Section 3 the 
infrared divergencies of the virtual loop amplitude are discussed. These di- 
vergences are cancelled either by real gluon emission or by the diagrams 
where an additional virtual gluon is emitted off the second t-channel (screen- 
ing) gluon in the whole pp — > p@bb®p amplitude. Numerical estimates are 
illustrated in Section 4. 

2 Calculation of the gg — > bb amplitude 

Here, and in what follows, we consider the amplitude, defined as 

T 99 ->bi = EE M^5 ab 5 ei62 , (9) 

a, b ei,<52 

where a and b are the colour indices (a, b = 1,2,...,N% — 1) and ei^ are 
the transverse polarization vectors of the incoming gluons. Introducing the 
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outgoing quark spinors, the amplitude can be written as 

Tgg^bb = u(Pl)T(^)v(p 2 ), (10) 

where the matrix T is built from the 4-momenta p involved in the reaction, 
f 1 = p^ . There are three independent 4-momenta, which we take to be pi, 
P2 of b and b, and the difference of the incoming gluon momenta, Pa~Pb- The 
matrix p'l can be moved to the left, and $ 2 to the right, until they disappear 
upon acting on the spinor, 

u(pi)j>\ = mu(pi), i> 2 v{p 2 ) = -mv(p 2 ), (11) 

where m = mt,. This rearrangement reduces the amplitude to two spinor 
structures, corresponding to the following helicity- violating and helicity-conserving 
terms^l 

T gg ^ b - b = u(p 1 )v(p 2 )T 1 + u(pi)(p ! a-P , b)v(p 2 )T 2 . (12) 
The scalar coefficients T 12 depend on the invariants 

s = (pi+p 2 ) 2 , t=(PA-Pi) 2 u={p A -p 2 ) 2 . (13) 

The function Ti is symmetric with respect to the t-u interchange, while T 2 is 
antisymmetric. Their calculation requires an evaluation of the spinor traces, 
which can easily be carried out with the help of an analytical program. 



2.1 The Born contribution 

As discussed above, only the helicity-violating piece contributes at the Born 

level to the J z = amplitude, that is the sum (++) + ( ) of the helicities 

of the incoming gluons, 



Born 



4nas Cp 2m 



m 2 — t 



+ 



77T 



T. 



Born 



0. 



(14) 



The loop integration is reduced to master scalar integrals, and the final re- 
sult is written in terms of the set of these integrals. The Feynman gauge 
is used in the MS-scheme with the dimensional parameter of order of the 



4 In principle, there are two other possible structures related to the fourth vector ~ 
s^ v \<jPi{pa -Pb) X P2 namely, u(pi)j 5 v(p 2 ) and u(pi)(i> A ~ i>B)l z v{p 2 ), but they do not 
appear for our subprocess. 
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characteristic virtualities of the process, 47r/z 2 = —k 2 = s/4. The ultraviolet 
(UV) divergencies are renormalized, to first order in as, by the gluon and 
quark Z-factors, mass and charge renormalization in the Born term. 

For central exclusive production, the effective luminosity of the incoming 
active gluons, 

M*dC(gg pp ) , 2 _ cl . . 

dydM 2 ~ b L ' (15) 

is given by the integral expression [2] 



m-\)h) q 



where the integration is over the transverse momentum of the gluon loop 
in Fig. [TJ Here b is the t-slope, corresponding to the momentum transfer 
distributions of the colliding protons, and x\p are the light-cone momentum 
fractions carried by the active gluons^]. The colour flow due to the active 
gluons is screened by the second t-channel gluon, which carries the momen- 
tum fractions x'f and x 2 ~ ■ The gap survival factor S 2 accounts for the soft 
rescattering effect, that is S 2 is the probability that the rapidity gaps are not 
populated by secondaries produced in possible soft rescattering, see [TJ [21 [21] • 
The integral over the gluon transverse momentum Q t is convergent, in both 
the ultraviolet and infrared regions. The conribution of the low Qt domain 
is suppressed by Sudakov-like form factors incorporated in the unintegrated 
gluon densities f g . 

Due to the presence of the second t-channel (screening) gluon, there is 
no infrared divergency in the matrix element of the 'hard' gg(J z = 0) — > bb 
subprocess. This opens up the possibility to regularize the infrared (IR) 
divergences by a fictitious gluon mass m g . In order to preserve gauge invari- 
ance, the mass m g can be formally introduced via the Higgs mechanisn0. 
Of course, this will generate additional diagrams with Higgs boson exchange 
or production. Such Higgs boson exchange will provide 'cross-talk' between 
the 'hard' matrix element and the second t-channel (screening) gluon. As a 
result, we set m g = Q t , that is, introduce a physical infrared cutoff in the 
exclusive bb production. After this, the contribution of the "artifical" Higgs 



The superscript + (— ) is to indicate the light-cone momentum fraction of the mo- 
mentum qi, ((72) of the first (second) colliding proton, that is Q M = x'^qx^ + £ 2 ~92/i +Qtfi- 
6 This way of introducing the infrared cutoff m g was used, for example, in the original 
BFKL calculations, see l25l. 
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boson diagrams can be eliminated by choosing the mass of the "Higgs" boson 
to be very large, M\ 3> s. Since there is no need to regularize the integral 
( I16p . here we will work in the D = 4 world, so the polarizations of our gluons 
are constrained to D = 4 space. 

2.2 One- loop effects 

Both the structures Ti j2 in ([12]) can contribute to the one-loop correction. 
We begin with the helicity-conserving part 

T 2 (s,t } u) = A 2 (s,t,u) - A 2 (s,u,t). 

We retain only the quark mass m = contribution. The function A 2 (s, t, u) 
is given by the set of Feynman diagrams shown in Fig. [2J The second term in 
the expression for T 2 (s, t, u) takes care of t <-> u interchange. When performed 
naively, a straightforward calculation results in a complete cancellation be- 
tween the two terms. So we have T 2 (s, t, u) = at the one- loop level, similar 
to the Born term0. 

The QHNC part of the amplitude is given by the two terms, 

Ti(s,*,u) = Ax(s,t,u) + Ai(s,«,t), (17) 

where the second term, as in the previous case, comes from the w-channel 
crossing. Recall that the amplitude T\ is symmetric with respect to the t-u 
interchange. As in the Born term, the one-loop amplitude vanishes if we set 
m = 0. Here we present an expression for the amplitude to lowest order in 
mass, that is 

A = 47ra 5 -Y-m [A + A q + A g + A qg ), (18) 

7 At first sight, it appears that the zero value of Ti is in contradiction with [26l I27j. 
which lead to a non-zero result for the one-loop (cut-non-reconstructible) QHC amplitude 
gg(Jz = 0) — > bb in dimensional regularization, assuming that the incomimg gluons are 
on-mass-shell. This non-zero cut-non-reconstructible contribution comes from a ratio of 
the form e/e = const, where the denominator (that is, the factor 1/e) is of infrared origin. 
However, in our case, due to the presence of the screening gluon, the (transverse) size of 
the interaction region is limited by a value of order of ~ 1/Qt- This generates a dynamical 
infrared cut-off at a scale ~ Q\ , which, in our calculation, is taken care of by introducing 
an effective gluon mass m g — Q t . As a result, the infrared singularities are absorbed in 
the unintegrated gluon structure functions. The factor 1/e is replaced by ln(/i/m ff ), and 
as e — * we obtain T2 = 0. Thus our procedure is quite different from the on-mass-shell 
calculations. 
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where m is the quark pole mass, m = Tn^rrib), and the running coupling 
constant is taken at the scale /x, as = «s(/i 2 ). The function A is made up of 
pieces coming from the logarithmically divergent scalar master integrals with 
two propagators. The 1/e terms, appearing in dimensional regularization, are 
subtracted by the counterterms coming from the field Z-factors, quark mass 
and coupling constant renormalization in the Born term. To be precise, 
A collects what is not included in the functions /, which result from the 
master integrals. It depends on which form has been assumed for the gluon 
polarizations in dimensional regularization. That is, whether the transverse 
tensor lies exactly in two-dimensional space D — 2 = 2, or whether it has 
an admixture of the extra dimension, D — 2 = 2 + 2e. As we discussed, in 
central diffractive production the active incoming gluons are screened (in the 
whole amplitude of the exclusive process) by a second t-channel gluon, see 
Fig. [H In this external loop we have neither infrared nor UV divergencies. 
Thus the polarisations of incoming gluons lie in two-dimensional D — 2 = 2 
space. Moreover, when using the unintegrated gluons obtained via the KMR 
prescription [2E], from the integrated MS gluon given by the global parton 
analyses, we have to account for the fact that the whole gluon renormalisation 
factor Zj, is already included in the incoming parton distribution. So the 
matrix element of the 'hard' subprocess must be calculated with the Z factors 
equal to 1 for all 'external' lines (just as is the case when the external lines are 
on- mass-shell) . In other words, here we present the result of the calculation^! 
with Z% = Z 2 = 1. 

The result obtained for Z2 = 1, Z% = 1, with the infinite parts subtracted, 
has the form, 



aD-2=2 _ ' 



-AC F N C In — + C F — c - In - — - - 2AC 2 F 



(19) 



The other functions in ([IB]) are given by the finite parts of the individual 
Feynman diagrams shown in Fig. [2] with the colour coefficients q, 



Cl 



c 2 



c 3 



r 2 



N 2 



2N r 



2N r 



1 N 2 - 

2N r 2N r 



8 The renormalisation factor, Z2, of the &-quark is included in the evolution which 
describes the fragmentation of the 6-quark jet. 
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c 5 


= iV c C F = 


c 6 


= N 2 C -1 


C 7 = Cg 




Cg = Cio 





N 2 c -l 

' 2Nr. 



2N C 
2 2N C 

The master integrals, ioooo( s >0 e ^ c ) appearing in the expressions below, are 
given in the Appendix A. We have 

A, = iCl , M s + Cn(t) 2[(2. + t)fe + C ,) + 2(2,- t)Cl -2e 4 ,] 

5 £ 



^9 = 4i oooo(M) c 5S + 3/f 000 (s)[c 6 - 2c 5 ] - 4c 5 / ? 100 (t) 

4[-(35 + 2t)c 5 + c 6 g] 4(2s + t)c 5 

"T J 1010l 6 J 2 J 0101V C J . 



^ -[2 S 2 + 3st + 3t 2 ](c 7 + c 9 ) A 3( S + t)(c 7 + c 9 ) 



s 



'/'/ 'iliiiiuS "• T ^ 1000 

, rgg m -3 (C7 + Cg) t (c 7 + Cg) (2s + 3t) 

+ J ioool r J : r ^ooioW : 

/ >\ (C7 + C 9 )(S + 3t) + 2(C 8 + Cip)g -4(C 7 + Cg) 

-'ooioW + i ioioW" 



+ 2/^10 (*) 



s s + t 

[(cs + c w )(t-2s) + 2(c 7 + c 9 )s] 



St 

3 The infrared contribution 

The amplitude Ti contains logarithmic infrared (IR) divergences which in 
the formulae presented above are regularized by the effective gluon mass m g 
cut-off. These IR divergent terms are essentially the usual Sudakov form 
factors, that is the probability not to emit additional gluons in the exclusive 
process pp — > p © bb © p. 

The Sudakov-like form factor, due to emission from the initial active 
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gluons, is equal to 

S{Q t ,n) = exp (- j* g|f zP gg {z)dz\ (20) 

with P gg {z) being the gluon-gluon Altarelli-Parisi LO splitting function cor- 
responding to real gluon emission, and A = k t /(fi + k t ). It was already 
included in the effective gluon-gluon luminosity C(gg pp ) used to calculate 
the exclusive cross section [2]. Therefore, we have to subtract the term 

/ r^ 2 rye rlh 2 rl-A 
nBorn / r a S al ^t 



T > [rkttl, zP " [z)dz ) (21) 

from the amplitude T. 

As discussed above, because of the presence of a screening gluon, the 
infrared cut-off is given by the transverse momentum (virtuality) of the in- 
coming active gluon, and in the expressions for the amplitudes A g and A qg , 
where the 'internal' gluon is radiated from a gluon line, we replace m g by Qt- 
This cancels the main (oc ln 2 (mp) part of the IR divergency 



T x ~ 2m(4vra s )iV c CV— [- + -) \n 2 (s/m 2 g ). (22) 

The logarithmic IR divergency in the amplitude A q is cancelled after account- 
ing for the real soft gluon emission in the 6-quark jet. Usually, in Monte Carlo 
simulations and/or jet searching algorithms, such emission is described by LO 
quark evolution. So we have to subtract the term 

^(-jCsfjT'H' (23) 

where now P qq {z) is the quark-quark splitting function; the lower limit k t $ 
is fixed by the experimental conditions - gluons with transerse momenta 
with respect to 6-jet axis k t < k t0 are included in the definition of the jet. 
The kinematic limit is A q = 2k t Q/y/s. As a result, we put m g = k t o in the 
expression^] for A q . 



9 There still remains a contribution proportional to the first power of ln(m g ), which 
is not cancelled by the subtractions (|21|23| . This contribution arises from large-angle 
soft-gluon emission, when we cannot neglect the interference between the emission from 
the gluon and from the quark lines. Such interference, hidden in the amplitude A qg , is 
not included, either in the definition of the jet or in the effective gg pp luminosity. The 
corresponding IR divergency is cut off by the presence of the screening gluon in the effective 
gg pp luminosity, that is by the gluon transverse momentum Q t ; to mimic this fact we set 
ra„ = Qt, as before. 
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4 The double-logarithmic contributions 



As was discussed in the Introduction, large double-logarithmic terms can 
be of Sudakov or non-Sudakov origin. The Sudakov contributions reflect 
the possibility to emit additional soft gluons. They are absorbed (and sub- 
tracted) in the definitions of the gg pp luminosity and in the prescription for 
the quark jet search. The non-Sudakov logarithms come from the kinemati- 
cal domain in which one of the quark propagators in the diagram is soft. In 
the case of the 77 — > bb process this contribution was numerically quite large 

T ncn-sud (77 ^ _ T? OTn ■ 3C F ^ In 2 (s/ m 2 b ), (24) 

4:71 

see (EJ). In our case, with a larger number of diagrams and larger colour 
coefficients in the one-loop gg — > bb Feynman graphs, there is a danger 
that we could find an even larger non-Sudakov DL correction. However, 
the situation appears to be different. Contributions which correspond to 
diagrams with 3 gluons in the loop and to diagrams with 2 gluons have 
different signs. This is analogous to the destructive interference between the 
emission of a photon from the incoming and the outgoing electron for small 
angle scattering. Unlike the 77 — > bb case, here we have additional double- 
logarithm contributions coming from A g , A gg , such that the final result does 
not contain a large numerical coefficient 

T™ n - Sud (gg bb) ~ T* orn ■ {3C F -N c )^ In 2 {s/m 2 b ). (25) 

47T 

We see that, instead of the naively expected factor 3N C (which indeed comes 
from the first term of A g , that is from the integral /oooo)) the coefficient in 
the sum of the A g and A qg amplitudes in front of the non-Sudakov double- 
logarithm is proportional to 3N C — 4N C = — N c . 



5 Discussion 



To evaluate the role of the one-loop correction numerically, we first calculate 
the cross section 

a (T (o)+(i) 1/1 _ \ 2 
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: Q=2 GeV, k t0 =2 GeV 



0.05 0.1 0.15 0.2 0.25 0.3 0.35 0.4 0.45 0.5 

cos9 

Figure 3: Angular dependence of the exclusive bb cross section for different 
choices of the infrared cutoffs. For comparison, we also show the Born result for 
two choices of the running 6-quark mass. 



that must be multiplied by the effective gluon-gluon luminosity (floT) [2 J . The 
results are shown in Fig. [3] for different values of the infrared cutoffs Q t and 
k t0 . The scale is taken to be Aufi 2 — s/A. 

The value Q t = 2 GeV corresponds to the maximum of the integrand in 
(fl6j) for the exclusive production of a Higgs boson of mass M H = 120 GeV at 
the LHC. The choice k t0 = 5 GeV appears to be reasonable for the standard 
6-quark-jet searching; the gluons with transverse momentum (with respect to 
6-jet axis) of k t > 5 GeV can be considered as separate jets. To demonstrate 
the dependence of the cross section 

a (o)+(i) on 

the values of Q t and k t o we 
also show the predictions for Q t = 6 GeV and k t0 = 2 GeV. We use a 6-quark 
pole mass m^ijn^) = 4.2 GeV and take as{Mz) = 0.118. For lower infrared 
cutoffs the probability not to emit an additional gluon decreases, and the 
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Figure 4: The scale dependence of the ratio of the NLO exclusive bb cross section 
to that calculated in Born approximation. 



cross section is smaller. 

For comparison, we also show in Fig. [3] cross sections calculated in the 
Born approximation with the same renormalisation scale (s/4) for the QCD 
as coupling and the 6-quark pole mass mi, taken at the same scale mbijfib) 
(upper curve). In addition we plot the Born result for a 6-quark mass taken 
at the scale y/s/2 (lower curve); this shows that a large part of the one- 
loop suppression of the cross section comes from the running of the 6-quark 
mass. However, note that for large-angle scattering we observe a stronger 
suppression of the cross section due to other radiative corrections. 

In Fig. H] we show the scale dependence of the ratio er(°) + M / cr(°) of the 
whole one- loop cross section, integrated over the region 60° < 9 < 120° (that 
is | cos 6*| < 1/2), to that calculated in Born approximation with the pole 
mass mb(mb) of 6-quark. The result is shown for two different masses of bb 
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Figure 5: The mass dependence of the ration of the NLO exclusive bb cross 
section to that calculated in Born approximation. 

system, namely M&& = \/s = 120 and 200 GeV. Here we put Q t = 2 GeV and 
k t0 = 5 GeV. It is seen that the scale dependence in the region of 4vT/i 2 ~ s/4 
is rather flat. 

Finally, in Fig. El we present the analogous ratio, a^ + ^ / a^°\ of the CED 
cross sections, integrated over the region of | cos6>| < 1/2, expected at the 
LHC for exclusive bb production in the central region (with the rapidity of the 
66-pair y = 0). In this case the gg pp — > bb amplitudes T® orn and T\, which 
enter the cross section ( 1261) . were convoluted with the luminosity amplitude 
(the integrand of ( |T6|) ) following the Q t factorisation prescription; that is the 
amplitude Ti(Q t ) was included inside the Q t integral in (fT6l) . Again we show 
the results for two values of the infrared cutoff in the 6-jet definition - k t o = 5 
GeV and k tQ = 2 GeV - as the function of the mass of the bb pair for the 
scale equal to M^/2. Of course, for the smaller value, k t0 = 2 GeV, we have 
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a stronger suppression, but this does not mean that by selecting narrower 
6-jets (with a smaller cone size AR or a smaller k t0 ) we can improve the 
Higgs signal-to-background ratio. The signal is diminished in the same way 
as the background when we suppress the emission of an additional gluon in 
the H — » bb decay; both the exclusive bb cross section and the H —> bb signal 
have the same Sudakov suppression. Thus, in order not to lose statistics, it 
is better not to take the value of k t0 to be too small. On the other hand, 
kto should not be too large. Otherwise, we will not sufficiently suppress the 
three-particle, gg pp — > bb + g, background |14j . From this viewpoint the 
predictions corresponding to k t o = 5 GeV look appropriate. 

In conclusion, the good news is that the radiative QCD (one-loop) correc- 
tions suppress the exclusive bb background (by a factor about 2, or more for 
larger 66-masses) for central exclusive diffract ive (CED) Higgs production, 
in comparison with that calculated using the Born gg pp — > bb amplitude. 
As discussed in [121 EJ; bb production in the \J Z \ = 2 state is another back- 
ground, which cannot be distinguished from the H — > bb decay. However, 
this contribution can be suppressed by selecting events with smaller trans- 
verse momenta of the forward outgoing protons [T2]. Therefore, the exclusive 
gluon-gluon dijet production becomes the most important background for the 
CED Higgs process, see [21 [6] for a detailed discussion. In order to suppress 
further this QCD background we need better experimental discrimination 
between 6-quark and gluon jets; that is, to achieve a lower probability P g /b 
for misidentifying a gluon as a fo-jet. 
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Appendix A: Master integrals 

We present below the scalar master integrals needed to compute the one-loop 
effects shown in Fig. [2j The results contain the functions 

m = f t dx \EH±*l = - ii2 (-«, 

Jo X 



p = \Js(s — Am 2 ), 

C = -- - In m2 

e 47r/i 2 ' 

where we work in D = 4 + 2e dimensions. 
Extra gluon between quark lines 

Here, we give the integrals for diagrams 1-4 of Fig. [2} which involve the 
virtualities 

a = k 2 —m 2 , b = {k—p B ) 2 —m 2 , c = (k+pA—pi) 2 —m 2 , d = {k+p A ) 2 —m 2 . 
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Extra gluon coupling to gluons 

Here, we give the integrals for diagrams 5 and 6 of Fig. [2j which involve the 
virtualities 

a 1 = k 2 -m 2 b' = (k-p B ) 2 -m 2 c' = (k+p A-Pif —m 2 , d! = (k+p A ) 2 -m 2 
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Extra gluon between a quark and a gluon line 

Here, we give the integrals for diagrams 7-10 of Fig. [21 which involve the 
virtualities 

a = k — to , b = (k+pb) —m, c = [k+p 1 —pA) —to , a = (k+pi) — m , 



^oooo( M 5 : _ / -9 ~r~7 
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Appendix B: MS renormalization 

It is useful to recall the well known one-loop renormalization in the MS- 
scheme. The bare quark mass to and coupling constant a® are expressed 
through their renormalized values to and as in the renormalized amplitude 
T R . In addition, the renormalization of the quark and gluon wavefunctions 
has to be taken into account, which results in factors Z 2 and Z3. The relation 
between the bare and physical coupling at one-loop order reads 



tt'o 



a 



1 + -(3 a + 0(a 2 ) 
e 



a + a, 



(1), 



a) 
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Similarly 



^2,3 



1 + Z 2,3> 



where Z^l are given by the first-order corrections to the quark and gluon 
propagators. The correction to the bare quark mass, 

m (m) = m + m^\m) = m + S(m), 

is determined by the first-order quark self energy S(^). With these factors 
included, the first-order amplitude takes the form 

t r = T onc-ioo P(ao5mo) + [ Z V + zi 1] +m£\m) + a£\a)] T Bom 



The quark one-loop self energy reads 

S(p) = -^ r ( 2 "f) (^[Z)m + x(2-^] 
where x = 1 — x and 7 is Euler constant. As a result 



—xxp 2 + xm 2 + xm 2 ' 
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The gluon polarization operator in Feynman gauge, 

rp» = (^V-p^)tt(p 2 ), 



reads to one-loop order 
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With D = 4 + 2e, one gets 
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